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Abstract. Using the new data from the OECD-WTO world network of economic activities we construct 
the Google matrix G of this directed network and perform its detailed analysis. The network contains 58 
countries and 37 activity sectors for years 1995, 2000, 2005, 2008, 2009. The construction of G, based on 
Markov chain transitions, treats all countries on equal democratic grounds while the contribution of activity 
sectors is proportional to their exchange monetary volume. The Google matrix analysis allows to obtain 
reliable ranking of countries and activity sectors and to determine the sensitivity of CheiRank-PageRank 
commercial balance of countries in respect to price variations and labor cost in various countries. We 
demonstrate that the developed approach takes into account multiplicity of network links with economy 
interactions between countries and activity sectors thus being more efficient compared to the usual export- 
import analysis. Our results highlight the striking increase of the influence of German economic activity 
on other countries during the period 1995 to 2009 while the influence of Eurozone decreases during the 
same period. We compare our results with the similar analysis of the world trade network from the UN 
GOMTRADE database. We argue that the knowledge of network structure allows to analyze the effects of 
economic influence and contagion propagation over the world economy. 

PACS. 89.75.Fb Structures and organization in complex systems - 89.65.Gh Econophysics - 89.75.He 
Networks and genealogical trees - 89.20.Hh World Wide Web, Internet 


1 Introduction 

The matrix analysis of Input-Out transactions had been 
pushed forward in the fundamental works of Leontief m 
becoming nowadays at the heart of modern treatment of 
economic relations (see e.g. [3]). The recent reports of 
the Organisation for Economic Co-operation and Devel¬ 
opment (OECD) [4] and of the World Trade Organization 
(WTO) [5] demonstrate the growing complexity of global 
manufactoring activities, exchange and trade in the mod¬ 
ern world. Thus the advanced matrix methods are highly 
desirable for the analysis of these complex systems. 

In parallel, during the last decade the modern society 
generated enormous communication and social networks 
including the World Wide Web (WWW), Wikipedia, Twit¬ 
ter and other directed networks (see e.g. Is])- The concept 
of Markov chains provides a powerful mathematical ap¬ 
proach for analysis of such networks. In particular, the 
PageRank algorithm, developed by Brin and Page in 1998 
[7] for the WWW information retrieval, became at the 
mathematical foundation of the Google search engine (see 
e.g. i). This algorithm constructs the Google matrix G 
of Markov chain transitions between network nodes and 
allows to rank billions of web pages of the WWW. The 
spectral and other properties of the Google matrix are 
analyzed in [9] . The history of the development of Google 


matrix methods and their links with research in social sci¬ 
ences and works of Leontief in economy is reviewed in in 

mm- 

The results presented in |12ll3j for the World Trade 
Network (WTN), constructed from the United Nations 
GOMTRADE database [14], show that the Google matrix 
analysis is well adapted to the ranking of world countries 
and trade products and to determination of the sensitivity 
of trade to price variations of various products. The new 
element of such an approach is an equal (“democratic”) 
treatment of all countries independently of their richness 
thus being rather different from the usual Import and Ex¬ 
port ranking. At the same time the contributions of var¬ 
ious products are considered being proportional to their 
trade value (volume) contribution in the exchange flows. 

In this work we use the Google matrix methods for 
analysis of the contagion effects on the World Network of 
Economic Activities (WNEA). We use the new database 
of the OECD-WTO with the network of 58 countries and 
37 activity sectors. At the difference of the World Trade 
matrix which report trade in goods between countries, 
the WNEA maps the imports and exports of intermediate 
goods and services between industries. Those globalised 
inter-industrial exchanges of intermediate inputs are one 
of the characteristics of the International Supply Chains 
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where the production of final goods results from the com¬ 
bination of various industrial tasks that are internationally 
outsourced. The first results of the Google matrix analy¬ 
sis of this database have been reported in |15ll6j for years 
1995 and 2008 while here we analyze the time evolution 
for years 1995, 2000, 2005, 2008, 2009. We show that our 
approach gives the results being different from the usual 
import-export flows for individual countries represented 


in Fig.[^ 
at dl]). 


for year 2005 (world map of countries is available 
The main reason of this difference is due to the 
fact that the Google matrix analysis takes into account the 
multiplicity of links between various nodes of the network 
while the import-export approach provides only local in¬ 
formation without taking into account the complex link 
relations between nodes. 

The new element of the OECD-WTO database is that 
it includes the transactions between different activity sec¬ 
tors while the COMTRADE database for multiproduct 
trade has no transitions between different products (even 
if they exist in reality, e.g. metal and plastic are used for 
production of cars). 

We note that there has been a number of other inves¬ 
tigations of the WTN reported in |18I19I2QI21I22I23I24] . 
However, in this work we have the new important ele¬ 
ments, developed in |12ll3ll5j : the analysis of PageRank 
and CheiRank probabilities corresponding to direct and 
inverted network flows and related to Import and Export; 
democratic treatment of countries combined with the con¬ 
tributions of sectors (or products) being proportional to 
their commercial exchange fractions. We point out that 
the OECD-WTO TiVA database of economic activities 
between world countries and activity sectors has been cre¬ 
ated very recently (2013) and thus this work represents 
new studies of the WNEA data evolving in time, extend¬ 
ing the results reported recently in ng. 


2 Methods and data description 

The list of Nc = 58 countries (57 plus 1 for the Rest Of 
the World (ROW) is given in Tablewith their flags. Fol¬ 
lowing [12] we use for countries ISO 3166-1 alpha-3 code 
available at Wikipedia. The list of sector activities with 
their names is given in Table . The sectors are classified 
according to the International Standard Industrial Classi¬ 
fication of All Economic Activities (ISIC) Rev.3 [2^ . 

For a given year, the TiVA data extend OECD In¬ 
put/Output tables of economic activity expressed in terms 
of USD for a given year. These data are tentative and had 
been released in 2013. A next version is expected to be 
available in 2015. From these data we construct the ma¬ 
trix Mcc' ss' of money transfer between nodes expressed 
in USD: ’ 



Fig. 1. World map of countries with color showing country im¬ 
port (top panel) and export (bottom panel) with economic ac¬ 
tivity (trade) value (volume) expressed in billions of USD and 
given by numbers at color bars; the gray color marks countries 
attributed to the ROW group (rest of the world) with exchange 
values 387 (Import) and 547 (Export) in billions of USD. The 
data are shown for year 2005 with Ac = 57 -h 1 countries (with 
ROW) for the economic activities in all Ng = 37 sectors. Coun¬ 
try names can be found in Table and in the world map of 
countries m- 


sector, a link gives a transfer from a sector of one country 
to another sector of another country. We construct the ma¬ 
trix Mcc',ss' from the OECD-WTO TiVA Input/Output 
tables using the transposed representation so that the vol¬ 
ume of products or sectors flows in a column from line to 
line. In the construction of Med ,ss' we exclude exchanges 
inside a given country in order to highlight the trade ex¬ 
change flows between countries (elements inside country 
are zeros). The method of construction of the Google ma¬ 
trix from the matrix Mcd,ss' are described in m (see m 
for the COMTRADE database) but for convenience of a 
reader we repeat this description here. 

We define the value of imports Ws and exports 17* for 
a given country c and sector s as 

^ , v: = V (2) 

d,s' d,s' 


Mcd,ss' = transfer from country c', sector s' to c, s (1) 

Here the country indexes are c, c' = 1,..., Ac and activ¬ 
ity sector indexes are s, s' = 1,..., A^ with Ac = 58 and 
Ng = 37. The whole matrix size is A = Ac x Ag = 2146. 
Here each node represents a pair of country and activity 


The import W = Zls ^cs and export values 

for countries c are shown on the world map of countries 
in Fig. I^for year 2005. We note that often one uses the 
notion of volume of export or import (see. e.g. [T3|) but 
from the economic view point it more correct to speak 
about value of export or import. 
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In order to compare later with the PageRank and 
CheiRank probabilities used below we define exchange 
value ranks in the whole matrix space of dimension N = 
NcX Ns. Thus the ImportRank (P) and ExportRank (P*) 
probabilities are given by the normalized import and ex¬ 
port values 


Pi = v,jv, p: = v:jv, ( 3 ) 

where i = s {c — l)A^s, i = and the total 

exchange value is V = ^c,c',s,s' ,ss' = ^c,s^cs = 

y R*. 

/-^CS cs 

The Google matrices G and G* are defined as N x N 
real matrices with non-negative elements: 


Gij — aSij^{l — a)viej , G'^ij — aS""ij-\-{l—a)v^ej , (4) 

where N = Nc x Ng, a G (0,1] is the damping fac¬ 
tor (0 < (a < 1), Cj is the row vector of unit elements 
{ej = 1), and Vi is a positive column vector called a per¬ 
sonalization vector with Vi = 1 |8ll3j . We note that 
the usual Google matrix corresponds to a personalization 
vector Vi = Ci/N with = 1. In this work, following 
PTTsl , we fix a = 0.5 noting that a variation of a in a 
range (0.5, 0.9) does not significantly affect the probability 
distributions of PageRank and CheiRank vectors msM- 
The choice of the personalization vector is specified below. 
Following m we call this approach the Google Personal¬ 
ized Vector Method (GPVM). 

The matrices S and S* are built from money matrices 

^^cc' ,ss' 


Si. 


s*,, 


1/N if V'.' = 0 
M,,c,s'slV:,,, if V;,, ^ 0 

1 /iv if yt,, = 0 


(5) 


where c, c' = 1,..., Nc] s, s' = 1,..., i = s-^{c—l)Ns] 
i' = s' (c' — l)Ns] and therefore i, = 1,..., A/". Here 
K's' = ^cc',ss'- The sum of elements of each column 

of S and S'* is normalized to unity and hence the matrices 
G,G*,S, S* belong to the class of Google matrices and 
Markov chains. Thus S, G look at the import perspective 
and S*, G* at the export side of transactions. 

PageRank and CheiRank (P and P*) are the right 
eigenvectors of G and G* matrices respectively at eigen¬ 
value A = 1. The equation for right eigenvectors have the 
form 


^ = • ( 6 ) 

3 3 

For the eigenstate at A = 1 we use the notation Pi = 
with the normalization ^ Pi = Xli = 1- 
For other eigenstates we use the normalization = 

Si = 1. The eigenvalues and eigenstates of G, G* 
are obtained by a direct numerical diagonalization using 
the standard numerical packages. 

The components of P^, P*^ are positive. In the WWW 
context they have a meaning of probabilities to find a ran¬ 
dom surfer on a given WWW node in the limit of large 


number of surfer jumps over network links [8]. In WNEA 
context nodes can be viewed and markets with a random 
trader transitions between them. We will use in the fol¬ 
lowing notation of network nodes. We define the PageR¬ 
ank K and CheiRank P* indexes ordering probabilities P 
and P* in a decreasing order as P(P) > P(P + 1) and 
P*(P) > P*(P* -h 1) with P,P* = 1,...,P. 

We note that the pair of PageRank and CheiRank vec¬ 
tors is very natural for economy and trade networks cor¬ 
responding to Import and Export flows. For the directed 
networks the statistical properties of the pair of such rank¬ 
ing vectors have been introduced and studied in |26l27ll2j . 

We compute the reduced PageRank and CheiRank prob¬ 
abilities of countries tracing probabilities over all sectors 
and getting Pc = Pcs = ^ («§ + (c - l)Ps) and 

P* = Zs ^cs = Zs ('^ + (^ “ with the corre¬ 

sponding Kc and P* indexes. In a similar way we obtain 
the reduced PageRank and CheiRank probabilities for sec¬ 
tors tracing over all countries and getting 
= Zc ^ + (c - ^)^s) = Zc Pcs and 

P* = Zc^* ('^ + (<^ “ = ^cPcs with their corre¬ 

sponding sector indexes Kg and P*. A similar procedure 
has been used for the multiproduct WTN data m- 

In summary we have Kg , P* = 1,..., Pg and Pc, P* = 

1,..., Pc- A similar definition of ranks from import and 
export exchange value can be done in a straightforward 
way via probabilities Pg , P*, Pc^ P* ? Pcs ? Pcs correspond¬ 
ing indexes P,, , Pc, P*, P, P* - 

To compute the PageRank and CheiRank probabili¬ 
ties from G and G*, keeping a “democratic”, or equal, 
treatment of countries (independently of their richness) 
and at the same time keeping the proportionality of ac¬ 
tivity sectors to their exchange value, we use the Google 
Personalized Vector Method (GPVM) developed in [13] 
with a personalized vector Vi in At the first iteration 
of Google matrix we take into account the relative prod¬ 
uct value per country using the following personalization 
vectors for G and G*: 




V* 




N^y, V*, ’ 

c /L^s' cs' 


(7) 


using the definitions (§ and the relation i = s-\-{c—l)Ng. 
This personalized vector depends both on sector and coun¬ 
try indexes. As for the multiproduct WTN in [13] we de¬ 
fine the second iteration vector being proportional to the 
reduced PageRank and CheiRank vectors in sectors, ob¬ 
tained from the GPVM Google matrix of the first itera¬ 
tion: 


v'ii) 




(8) 


In this way we keep democracy in countries but keep con¬ 
tribution of sectors proportional to their exchange value. 
This second iteration personalized vectors are used in the 
following computations and operations with G and G* giv¬ 
ing us the PageRank and CheiRank vectors. This proce¬ 
dure with two iterations forms our GPVM approach. The 
difference between results obtained from the first and sec¬ 
ond iterations is not very large, but the personalized vector 
for the second iteration gives a reduction of fluctuations 











4 


V.Kandiah, H.Escaith and D.L.Shepelyansky: Google matrix of the world network of economic activities 


m- Below, in all Figures we show the GPVM results after 
the second iteration. 

As for the WTN it is convenient to analyze the distri¬ 
bution of nodes on the PageRank-CheiRank plane A"*). 
In addition to two ranking indexes K, if* we use also 
2DRank index K 2 which describes the combined contri¬ 
bution of two ranks as described in [27]. The ranking list 
K 2 {i) is constructed by increasing K ^ K+1 and increas¬ 
ing 2DRank index K 2 {i) by one if a new entry is present 
in the list of first if* < K entries of CheiRank, then the 
one unit step is done in iC* and K 2 is increased by one if 
the new entry is present in the list of first K < iC* entries 
of CheiRank. More formally, 2DRank 7^2 (^) gives the or¬ 
dering of the sequence of nodes, that appear inside the 
squares [1,1; K = = k; ...] when one runs progres¬ 

sively from k = 1 to N. Additionally, we analyze the dis¬ 
tribution of nodes for reduced indexes {Kc^ 


3 Results 

We apply the GPVM approach to the data sets of OECD- 
WTO TiVA of WNEA and present the obtained results 
below. 


Sectors (GPVM, 2005) 


Sectors (GPVM,2009) 
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Fig. 2. Gountry positions on PageRank-GheiRank plane 
(Kc,Kc) obtained for the WNEA by the GPVM analysis (top 
panels), ImportRank-ExportRank of exchange value/volume 
(bottom panels). Left (right) panels show year 2005 (2009). 


PageRank-CheiRank plane (Kc^K*) presented in Fig. 
for WNEA in years 2005, 2009. In the same figure we 
present the rank distributions obtained from ImportRank- 
ExportRank probabilities of exchange value. For the GPVM 
data we see the global features already discussed in |12ll5j : 
the countries are distributed in a vicinity of diagonal Kc = 
since for each country the size of imports is correlated 
with the size of exports, even if trade is never exactly 
balanced and some countries can sustain significant trade 
surplus or deficit. The top 20 list of top K 2 countries re¬ 
cover 13 of 19 countries of G20 major world economies (EU 
is the number 20) thus obtaining 58% (2005) and 63% 
(2009) of the whole list. This is close to the percentage 
68% obtained in na for year 2008. The Google ranking 
for WNEA (top panels in Fig. gives different positions 
for specific countries. Thus Russia and Saudi Arabia go 
to top Kc index values in PageRank comparing to Impor- 
tRank showing that their economies are highly sensitive 
to activity of petroleum sector. Similar features for these 
two countries are visible in 1995, 2008 m- 


PageRank-CheiRank (GPVM) ImportRank-ExportRank 
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• 1C01T05AGR 

• 2C10T14MIN 

• 3C15T16FOD 

• 4C17T19TEX 

• 5C20WOD 

• 6C21T22PAP 
7 C23 PET 

■ 8C24CHM 

■ 9C25RBP 

■ 10 C26 NMM 


■ 11C27MET 

■ 12 C28 EBM 

■ 13 C29 MEQ 
14 C30 ITQ 

▲ 15C31ELQ 

▲ 16C32CMQ 
A 17C33SCQ 
A 18C34MTR 

▲ 19C35TRQ 

▲ 20C36T37OTM 


21C40T41EGW ^ 31C73RDS 
-1- 22C45CON ^ 32C74BZS 
-1- 23C50T52WRT ^ 33 C75 GOV 
4- 24C55HTR ^ 34C80EDU 
-1- 25 C60T63 TRN 35 C85 HTH 
-1- 26C64PTL ♦ 36C90T93OTS 
-h 27C65T67FIN ♦ 37 C95 PVH 

28 C70 REA 
^ 29C71RMQ 
^ 30C72ITS 


3.1 Ranking of countries and sectors 

After tracing the probabilities P{K), {K"^) over sec¬ 

tors we obtain the distribution of world countries on the 


Fig. 3. Two-dimensional ranking of sectors on the (KsKs) 
plane using the GPVM approach for PageRank and GheiRank 
(left panels) and ImportRank-ExportRank (right panels). Each 
sector is represented by its specihc combination of color and 
symbol. The list of all 37 sectors are given in Table Top 
panels are for year 2005 and bottom panels are for year 2009. 
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After tracing over countries we obtain the PageRank- 
CheiRank plane {Kg, Kg* of activity sectors shown in Fig.^ 
As in na we find that some sectors are export oriented 
(e.g. s = 2 CIOTI 4 Mining at = 1 in 2009) others 
are import oriented (e.g. s = 22 C45 CON Construc¬ 
tion at Kg = 1 in 2009). The ImportRanking gives a 
rather different import leader s = 7 C23 Manufacture of 
coke, refined petroleum products etc. with Kg = 1 in 2009. 
Thus the Google ranking highlights highly connected net¬ 
work nodes while Import-Export gives preference to high 
value neglecting existing network relations between vari¬ 
ous countries and activity sectors. 

3.2 Price shocks and trade balance sensitivity 

On the basis of the obtained WNEA Google matrix we can 
now analyze the trade balance in various activity sectors 
for all world countries. Usually economists consider the 
export and import of a given country as it is shown in 
Fig. a Then the trade balance of a given country c can be 
defined making summation over all sectors: 

Be = J2iPes-Pcs)/J2(P:s+Pcs) = {P^ - Pe) / + Pc) ■ 

s s 

(9) 

In economy, Pc^Pf are defined via the probabilities of 
trade value Pcsj Pfg from (^. In our matrix approach, we 
define Pcs^Pfg as PageRank and CheiRank probabilities. 
In contrast to the Import-Export value our approach takes 
into account the multiple network links between nodes. 

The comparison of CheiRank-PageRank balance with 
Export-Import balance for the world countries is shown 
in Fig. for year 2009. Each country is shown by color 
which is proportional to the country balance Be (§ with 
the color bar given on the figure. For Export-Import bal¬ 
ance we see the dominance of petroleum producing coun¬ 
tries Saudi Arabia, Russia, Norway with the largest values 
in 2009 and 2008 (see Fig. 13 in [l5]). The CheiRank- 
PageRank balance highlights new features placing on the 
top Russia, Norway, Germany, China in 2008 m- In fact 
in 2008, USA has now a slightly positive balance in top 
panel of Fig. 13 while it was negative before in bottom 
panel of same figure. In 2009 after the world crisis there 
is a significant change for CheiRank-PageRank balance in 
the top panel of Fig.|^in 2009: USA takes the leading posi¬ 
tion while Saudi Arabia becomes even negative. The vari¬ 
ation of CheiRank-PageRank balance ABc = 5c(2009) — 
5c(2008) from 2008 to 2009 is shown in Fig.^ The strongest 
positive variation is obtained by Ireland, U5A and Japan, 
the strongest negative variation is for Saudi Arabia and 
Norway. We see that the broad network of economic ac¬ 
tivity relations and links makes the economies of the above 
countries more important in the world economy while Saudi 
Arabia, with the largest positive Export-Import balance, 
looses its leading position. 

To analyze the sensitivity of price variation in a cer¬ 
tain activity sector s we increase from 1 to 1 + the 
money transfer in the sector s in Mecss' in Q, where 5g 
is a dimensionless fraction variation of price in this sector. 



Fig. 4. World map of CheiRank-PageRank balance 5c = 
(PJ ~ Pc )/(PJ + Pc) determined for all Nc = 58 countries in 
year 2009. Top panel shows the probabilities P and P* given 
by PageRank and CheiRank vectors; the value of ROW group 
is 5 c= 58 = —0.02 02. Bottom panel shows the probabilities P 
and P* computed from the Export and Import value; the value 
of ROW group is 5c=58 = 0.0637. Names of the countries are 
given in Table and in the world map of countries m- 



Fig. 5. Difference of CheiRank-PageRank balance ABc = 
5c(2009) — 5c(2008) between years 2009 and 2008 shown by 
color for the world countries; for the ROW group we have 
ABc = —0.043 (gray). Names of the countries can be found 
in Tableland in the world map of countries Ha- 


After that the matrices G, G* are recomputed in the usual 
way described above and their rank probabilities P, P* are 
determined. Then we compute the derivatives of probabil¬ 
ities balance dBejddg over a price variation dg in a specific 
sector s. Of course, the computation is done at small val¬ 
ues of dg when the derivative is independent of dg and all 
price variations are kept sufficiently small. 
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Fig. 6. Derivative of probabilities balance dBc/dSr over price 
of sector s — 7 C23PET for year 2009. Top panel shows the 
case when Be is determined by CheiRank and PageRank vec¬ 
tors as in the top panel of Fig|3 the value of ROW group is 
dB^s/dSr = 0.0414. Bottom panel shows the case when Be 
is computed from the Export-Import value as in the bottom 
panel of Eig[^ the value of ROW group is dB^g/dSr = —0.0637. 
Names of the countries can be found in Table[T]and in the world 
map of countries CU- 



Fig. 7. Derivative of probabilities balance dBejdae' over la¬ 
bor cost of China c = 37 for year 2009. Top panel shows 
the case when Be is determined by CheiRank and PageRank 
vectors; here the special values are dB^s/dasr = —0.0163 for 
ROW group (gray) and dB^r/da^r = 0.3253 for China (ma¬ 
genta) . Bottom panel shows the case when Be is computed from 
the Export-Import value; the special values are dB^^jda^i — 
—0.0381 for ROW group (gray) and dRsr/dcrsr = 0.4732 for 
China (magenta). Names of the countries can be found in Table 
[^and in the world map of countries [17j . 


The sensitivity of country balance dBc/ddj to price 
variation of sector s = 7 Manufacture of coke, refined 
petroleum products and nuclear fuel is shown in Fig. 

For Export-Import in bottom panel the most sensitive 
countries are Lithuania (positive) and Vietnam (negative). 
Lithuania does not produce petroleum, but in fact in 2008 
there was a large oil refinery company there which had a 
large exportation value (see e.g. 

http://en.wikipedia.org/wiki/Economy_of_Lithuania). The 
Export-Import approach shows that Russia is slightly pos¬ 
itive, even less positive is Saudi Arabia, China and Ger¬ 
many are close to zero change, USA is only very slightly 
positive. The results of CheiRank-PageRank sensitivity 
(top panel) are significantly different showing strongly pos¬ 
itive sensitivity for Saudi Arabia, Russia and strongly neg¬ 
ative sensitivity for China, Germany and Japan; USA goes 
from slightly positive side in bottom panel to moderate 
negative one in top panel. The CheiRank-PageRank bal¬ 
ance demonstrates much higher sensitivity of Russia, Saudi 
Arabia and China to price variations of 5 = 7 sector com¬ 
paring to the case of Export-Import value analysis. The 
economies of Germany, China and Japan are also very 
sensitive to petroleum prices that is correctly captured by 
our analysis. We consider that the CheiRank-PageRank 
approach describes the economic reality from a new com¬ 


plementary angle and that provides new useful informa¬ 
tion about complex trade systems. We also note that the 
highly negative sensitivity of China to petroleum prices 
has been also obtained on the basis of Google matrix anal¬ 
ysis of COMTRADE data (see Fig.21 in pR]). 

It is also possible to determine the partial balance for 
a given sector s and given country c and to study its sensi¬ 
tivity to price variations in a sector s'. We do not discuss 
these characteristics here but an interested reader can find 
these results for year 2008 in m- 

Of course, the above derivatives over price of activity 
sector and labor country cost give only an approximate 
consideration of effects of price variations which is a very 
complex phenomenon. For an economic discussion of the 
effect of price shocks on international production networks 
we address a reader to the research performed in [28] . We 
will see below that our approach gives results being in 
a good agreement with economic realities thus opening 
complementary possibilities of economic activity analysis 
based on the underlying network relations between coun¬ 
tries and activity sectors which are absent in the usual 
Import-Export consideration. We present the results on 
sensitivity to sector prices and labor cost in next subsec¬ 
tions. 
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Fig. 8. Time evolution of the derivative dBajdac' over the 
labor cost c = 11 of Germany for years 1995, 2000, 2005, 
2008, 2009. For these years the special values are respectively: 
dB^sIdaii = -0.0402, -0.0307, -0.0351, -0.0367, -0.0388 
for ROW group (gray); dBwjdaw — 0.33, 0.3274, 0.3290, 
0.3248, 0.3760 for Germany (magenta). Names of the countries 
can be found in Table and in the world map of countries CD- 



Fig. 9. Time evolution of the derivative dBcjdac' over the la¬ 
bor cost of Eurozone monetary union (state in 2008 composed 
of 15 countries) for years 1995, 2000, 2005, 2008, 2009. For 
these years the special values are respectively: dB^g/daez = 
-0.1259, -0.0957, -0.0992, -0.0908, -0.0921 for ROW group 
(gray); and dBez/daez = 1.8422, 1.9235, 1.9394, 1.9418, 1.9508 
for Eurozone (magenta). Eurozone is composed from 15 coun¬ 
tries of its state in 2008: Austria (c = 2), Belgium (3), Finland 
(9), France (10), Germany (11), Greece (12), Ireland (15), Italy 
(17), Luxembourg (20), Netherlands (22), Portugal (26), Slove¬ 
nia (28), Spain (29), Malta (56). Names of the countries can 
be found in Table and in the world map of countries CD- 
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3.3 World map of sensitivity to labor cost 

The new element of the WNEA, compared to the multi¬ 
product WTN, is existence of transfers between sectors of 
the same economy. This allows us to consider the sensitiv¬ 
ity not only to sectoral prices but also the sensitivity to 
labor cost in a given country c (e.g. price shock affecting 
all industries in the same country). This can be taken into 
account by the introduction of the dimensionless labor 
cost change in a given country c by replacing the related 
monetary flows from coefficient 1 to 1 + dc in Mcc',ss] Q 
for a selected country c. 

Using the established structure of WNEA we can study 
the sensitivity of country balance dBc/dcr'^ to the labor 
cost in different countries. At the difference of sectoral 
shocks on one product, here the price shock affects all in¬ 
dustries in a country. As before, the change in price has to 
be small enough for the resulting simulation to remain in 
a neighbourhood of the original data. In fact we compute 
numerically the derivate dBcjda'^ corresponding only to 
small values of cr^. Indeed, larger shocks would trigger a 
series of substitution effects diverting trade to other part¬ 
ners. The modelling in the case of large shock variations 
is a very difficult task (see discussion and analysis of such 
situations at [29]). 

The derivative dBcjda'^ is computed numerically as 
described above. The world sensitivity to the labor cost 
of China is shown for year 2009 in Fig. Of course, the 
largest derivative is found for China itself {dBc/dac at 
c = 37 from Table [^. The effect on other countries is 
given by non-diagonal derivatives at c 7 ^ c' = 37. From 
the CheiRank-PageRank balance we find that the most 
strong negative effect (minimal negative dBc/dac') is ob¬ 
tained for USA, Germany, UK; a positive derivative is visi¬ 
ble only for Chinese Taipei {s = 38) and S.Korea {s = 19). 
For the Export-Import balance the results are rather dif¬ 
ferent: at first all derivatives at c 7 ^ c' are negative; among 
the most negative values are such countries as Hong Kong 
(most negative with dark red color but hardly visible due 
to its small size), Chinese Taipei, S.Korea, Vietnam. This 
sensitivity map for year 2008 is given in m, it has rather 
similar features. Thus the Google matrix approach brings 
a new perspective for analysis of complex of economic re¬ 
lations between countries and sectors. 

Another results for the time evolution of effects of la¬ 
bor cost in Germany are shown in Fig. [^for all available 
years from 1995 to 2009. Here we present results only for 
CheiRank-PageRank balance since the results of Export- 
Import balance, presented for year 2008 in m do not 
capture efficiently the multiple link relations. This time 
evolution demonstrates a spectacular increase of German 
influence from 1995 to 2009. Indeed, in 1995 the most 
strong negative sensitivity to German labor cost is visible 
mainly for USA but with time the influence of German 
economic activity extends to Russia and China captur¬ 
ing the large fraction of the whole world. Inside EU this 
influence is maximally negative in 2005 but it is slightly 
reduced in 2008 and 2009. 

For comparison we present the same time evolution 
of the sensitivity to the labor cost in Eurozone in Fig. 


Here, Eurozone is composed by 15 countries present in the 
Euro monetary union in 2008 and the labor variation ac' 
is taken to be the same in all these countries when com¬ 
puting the derivative dBc/dac' for all available years. It is 
striking to see that the evolution of Eurozone influence is 
opposite to Germany. Thus in 1995 we see a strong influ¬ 
ence on USA which however decreases significantly in 2008 
and 2009. The Eurozone influence on Russia (even if not 
so strong as for USA) also decreases with time. There is 
only a certain influence increase on China which however 
remains steady for years 2005 - 2009. 

4 Discussion 

In this work we have developed the Google matrix anal¬ 
ysis of the world network of economic activities from the 
OECD-WTO TiVA database. The PageRank and CheiRank 
probabilities allowed to obtain ranking of world countries 
independently of their richness being mainly determined 
by the efficiency of their economic relations. The devel¬ 
oped approach demonstrated the asymmetry in the eco¬ 
nomic activity sectors some of which are export oriented 
and others are import oriented. 

The CheiRank-PageRank balance Be allows to deter¬ 
mine economically rising countries with robust network 
of economic relations. The sensitivity of this Be to price 
variations and labor cost in various countries determines 
the hidden relations between world economies being not 
visible via usual Export-Import exchange analysis. 

Our Google matrix analysis highlights the striking in¬ 
crease of the influence of German economic activity on 
the economy of world countries during the period 1995 to 
2009. At the same time the influence of Eurozone decreases 
significantly. 

The knowledge of network connections in WNEA al¬ 
lows to investigate contagion propagation over the whole 
world. Indeed, a significant increase of petroleum prices 
can produce a shock wave which will propagate over the 
most sensitive links highlighted in our studies. We note 
that we consider the case of price contagion effect. We 
plan to develop and study the models of such shock con¬ 
tagion propagation in the future works. 

The comparison with the multiproduct world trade 
network from UN COMTRADE shows certain similari¬ 
ties between the two networks of WNEA and WTN. At 
the same time the WNEA data provides new elements 
for interactions of activity sectors while there are no di¬ 
rect interactions of products in COMTRADE database. 
From this viewpoint the OECD-WTO data captures the 
economic reality on a deeper level. But at the same time 
the OECD-WTO network is less developed compared to 
COMTRADE (less countries, years, sectors). Thus it is 
highly desirable to extend the OECD-WTO database. 

We think that the Google matrix analysis developed 
here and in fT 2 lT^ captures better the new reality of mul¬ 
tifunctional directed tensor interactions and that the uni¬ 
versal features of this approach can be also extended to 
multifunctional financial network flows which now attract 
an active interest of researchers . Unfortunately, the 
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data on financial flows have much less accessibility com¬ 
pared to the networks discussed here. 
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country name 

country code 

country flag 


country name 

country code 

country flag 

1 

Australia 

AUS 


30 

Sweden 

SWE 


2 

Austria 

AUT 


31 

Switzerland 

CHE 

o 

3 

Belgium 

BEL 

U 

32 

Turkey 

TUR 


4 

Canada 

CAN 


33 

United Kingdom 

GBR 


5 

Chile 

CHL 


34 

United States 

USA 

# 

6 

Czech Republic 

CZE 


35 

Argentina 

ARG 

/< 

o ^ 

7 

Denmark 

DNK 


36 

Brazil 

BRA 


8 

Estonia 

EST 


37 

China 

GHN 


9 

Finland 

FIN 


38 

Chinese Taipei 

TWN 

• 

10 

France 

FRA 

o 

39 

India 

IND 

ASk. 

11 

Germany 

DEU 

• 

40 

Indonesia 

IDN 


12 

Greece 

GRC 


41 

Russia 

RUS 

IP 









13 

Hungary 

HUN 


42 

Singapore 

SGP 


14 

Iceland 

ISL 


43 

South Africa 

ZAP 


15 

Ireland 

IRL 

o 

44 

Hong Kong 

HKG 

Q 

16 

Israel 

ISR 

✓ X 

45 

Malaysia 

MYS 


17 

Italy 

ITA 

o 

46 

Phillippines 

PHL 

d 

18 

Japan 

JPN 


47 

Thailand 

THA 


19 

Korea 

KOR 


48 

Romania 

ROU 

U 

20 

Luxembourg 

LUX 

s 

49 

Vietnam 

VNM 

o 

21 

Mexico 

MEX 

u 

50 

Saudi Arabia 

SAU 

Q 









22 

Netherlands 

NLD 


51 

Brunei Darussalam 

BRN 


23 

New Zealand 

NZL 

• 

52 

Bulgaria 

BGR 

w 

24 

Norway 

NOR 


53 

Cyprus 

CYP 










25 

Poland 

POL 


54 

Latvia 

LVA 


26 

Portugal 

PRT 


55 

Lithuania 

LTU 

• 

27 

Slovak Republic 

SVK 


56 

Malta 

MLT 


28 

Slovenia 

SVN 


57 

Cambodia 

KHM 

Q 

29 

Spain 

ESP 

e 

58 

Rest of the World 

ROW 

9 


Table 1. List of Nc = 58 countries (with rest of the world ROW) with country name, code and flag. 
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OECD ICIO Category 

ISIC Rev. 3 correspondence 

1 

C01T05 AGR 

01 - Agriculture, hunting and related service activities 

02 - Forestry, logging and related service activities 

05 - Fishing, operation of fish hatcheries and fish farms; service activities incidental to fishing 

2 

C10T14 MIN 

10 - Mining of coal and lignite; extraction of peat 

11 - Extraction of crude petroleum and natural gas; service activities incidental to oil and gas extraction excluding surveying 

12 - Mining of uranium and thorium ores 

13 - Mining of metal ores 

14 - Other mining and quarrying 

3 

C15T16 FOD 

15 - Manufacture of food products and beverages 

16 - Manufacture of tobacco products 

4 

C17T19 TEX 

17 - Manufacture of textiles 

18 - Manufacture of wearing apparel; dressing and dyeing of fur 

19 - Tanning and dressing of leather; manufacture of luggage, handbags, saddlery, harness and footwear 

5 

C20 WOD 

20 - Manufacture of wood and of products of wood and cork, except furniture; 

Manufacture of articles of straw and plaiting materials 

6 

C21T22 PAP 

21 - Manufacture of paper and paper products 

22 - Publishing, printing and reproduction of recorded media 

7 

C23 PET 

23 - Manufacture of coke, refined petroleum products and nuclear fuel 

8 

C24 CHM 

24 - Manufacture of chemicals and chemical products 

9 

C25 RBP 

25 - Manufacture of rubber and plastics products 

10 

C26 NMM 

26 - Manufacture of other non-metallic mineral products 

11 

C27 MET 

27 - Manufacture of basic metals 

12 

C28 FBM 

28 - Manufacture of fabricated metal products, except machinery and equipment 

13 

C29 MEQ 

29 - Manufacture of machinery and equipment n.e.c. 

14 

C30 ITQ 

30 - Manufacture of office, accounting and computing machinery 

15 

C31 ELQ 

31 - Manufacture of electrical machinery and apparatus n.e.c. 

16 

C32 CMQ 

32 - Manufacture of radio, television and communication equipment and apparatus 

17 

C33 SCQ 

33 - Manufacture of medical, precision and optical instruments, watches and clocks 

18 

C34 MTR 

34 - Manufacture of motor vehicles, trailers and semi-trailers 

19 

C35 TRQ 

35 - Manufacture of other transport equipment 

20 

C36T37 OTM 

36 - Manufacture of furniture; manufacturing n.e.c. 

37 - Recycling 

21 

C40T41 EGW 

40 - Electricity, gas, steam and hot water supply 

41 - Collection, purification and distribution of water 

22 

C45 CON 

45 - Construction 

23 

C50T52 WRT 

50 - Sale, maintenance and repair of motor vehicles and motorcycles; retail sale of automotive fuel 

51 - Wholesale trade and commission trade, except of motor vehicles and motorcycles 

52 - Retail trade, except of motor vehicles and motorcycles; repair of personal and household goods 

24 

C55 HTR 

55 - Hotels and restaurants 

25 

C60T63 TRN 

60 - Land transport; transport via pipelines 

61 - Water transport 

62 - Air transport 

63 - Supporting and auxiliary transport activities; activities of travel agencies 

26 

C64 PTL 

64 - Post and telecommunications 

27 

C65T67 FIN 

65 - Financial intermediation, except insurance and pension funding 

66 - Insurance and pension funding, except compulsory social security 

67 - Activities auxiliary to financial intermediation 

28 

C70 REA 

70 - Real estate activities 

29 

C71 RMQ 

71 - Renting of machinery and equipment without operator and of personal and household goods 

30 

C72 ITS 

72 - Computer and related activities 

31 

C73 RDS 

73 - Research and development 

32 

C74 BZS 

74 - Other business activities 

33 

C75 GOV 

75 - Public administration and defense; compulsory social security 

34 

C80 EDU 

80 - Education 

35 

C85 HTH 

85 - Health and social work 

36 

C90T93 OTS 

90 - Sewage and refuse disposal, sanitation and similar activities 

91 - Activities of membership organizations n.e.c. 

92 - Recreational, cultural and sporting activities 

93 - Other service activities 

37 

C95 PVH 

95 - Private households with employed persons 


Table 2. List of sectors considered by Input/Output matrices from OECD database, their correspondence to the ISIC classifi¬ 
cation is also given. 




